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THEORETICAL  CALCULATIOIIS  ON  THE  SHOCK  INITIATION  OP 

LIQUID  TNT 


By 

JULIUS  V;.  ENIG 
FRED  T.  nSTCALF 


ABSTRACT:  Theoretical  calculations  describing  the  initiation 
phenomena  in  honogeneons  liquid  TITT,  produced  by  shocks 
between  80  and  89  kilobars,  are  given.  The  initiation  mechan¬ 
ism  and  grovith  to  detonation  are  shovm  to  be  in  semi-quantita¬ 
tive  agreement  with  experimental  results.  Upon  entering  the 
explosive,  the  shock  wave  initiates  chemical  reaction.  After 
an  induction  time,  the  explosive  cell  that  has  been  heated  the 
longest  detonates.  The  detonation  vjave,  traveling  in  pre¬ 
compressed  explosive  v;ith  a  velocity  greater  than  the  steady 
state  velocity  in  uncompressed  explosive,  overtakes  the  decay¬ 
ing  initial  shock  and  temporarily  overdrives  detonation  in  the 
uncompressed  explosive.  The  numerical  experiments,  which  are 
based  on  finite  difference  solutions  of  the  partial  differen¬ 
tial  equations  of  hydrodynamics,  show  the  extreme  sensitivity 
of  the  induction  time  to  the  initiating  pressure  in  idealized 
one-dimensional  gap  tests.  Of  particular  Interest  is  the 
situation  at  onset  of  detonation  after  critical  initiation  by 
the  attenuated  shock  from  a  short  donor.  In  this  case,  the 
Interface  pressure  ir.  lovjer  than  that  in  the  shock  front 
although  the  Interface  temperature  is,  of  course,  higher  than 
that  of  the  shock  front. 
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INTRODUCTION 

The  initiation  of  detonation  in  explosives  by  shocks  has 
been  the  subject  of  a  very  large  number  of  experimental  investi¬ 
gations  and  a  growing  number  of  theoretical  calculations.  These 

1  2 

are  well  summarized  in  reviews  by  Bowden  and  Yoffee  ,  Jacobs  , 

^  2i 

Evans  and  Ablow"^,  and  Macek  .  Shock  Initiation  in  liquid  and 
solid  explosives  can  be  described  by  physico-chemical  inter¬ 
actions  that  occur  in  homogeneous  and  heterogeneous  materials, 
respectively.  Therefore,  a  theoretical  description,  while 
difficult,  is  still  simpler  for  homogeneous  liquids  than  for 
heterogeneous  solids  and  will  be  the  purpose  of  this  paper.  In 
the  following  paper^,  theoretical  calculations  for  shock  initia¬ 
tion  of  detonation  in  heterogeneous  solid  explosives  are  given. 

Hubbard  and  Johnson^  studied  the  properties  of  the  shock 
necessary  to  initiate  detonation  in  homogeneous  explosives  and 
found  that  when  an  element  of  explosive  was  subjected  to  a 
shock  wave  of  proper  strength  a  time  Interval  of  shock  heating 
existed,  called  the  "induction  time",  in  which  very  little  chemi¬ 
cal  reaction  occurred.  Then  complete  reaction  occurred  in  a 
relatively  short  time  and  a  wave  of  chemical  reaction  moved 
forward  as  a  detonation  wave.  In  particular,  they  gave  numeri¬ 
cal  solutions  to  the  one -dimensional  equations  representing  the 
conservation  of  mass,  momentum,  and  energy,  the  equations  of 
state,  and  the  equation  of  first-order  chemical  kinetics  for  a 
rectangular  pressure  pulse  applied  to  the  surface  of  a  semi- 
infinite  slab.  Prom  these  results  they  drew  conclusions  as  to 
the  relationship  between  the  magnitude  of  the  shock  and  the 
duration  of  the  pulse  that  led  to  initiation.  The  method  of 
solution  Involved  the  use  of  the  "q"  method  of  von  Neumann  and 
Rlchtmyer^.  They  assumed  the  same  Abel  tsTpe  equation  of  state 
to  hold  for  both  unreacted  explosive  and  product  gas . 
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Enig  obtained  numerical  solutions  of  the  one -dimensional 
Navier-Stokes  equations,  including  a  first-order  chemical 
reaction  term  of  the  Arrhenius  type,  for  the  problem  of  a  piston 
pushing  with  constant  velocity  against  one  surface  of  a  chemi¬ 
cally  inert  slab  whose  other  surface  is  in  intimate  contact  with 
a  semi-infinite  homogeneous  explosive.  The  inert  and  the  explo¬ 
sive  were  chosen  to  have  the  same  equation  of  state.  If  the 
shock  which  is  transmitted  into  the  explosive  from  the  inert  is 
sufficiently  strong,  chemical  reaction  is  initiated  at  the 
interface  between  the  inert  and  the  explosive.  By  choosing  the 
appropriate  piston  velocity  it  was  possible  to  produce  solutions 
in  which  the  shock  had  been  transmitted  some  distance  into  the 
explosive  before  the  inert-explosive  interface,  which  had  been 
heated  and  therefore  reacting  the  longest,  detonated.  This 
detonation  wave,  traveling  over  shocked  explosive,  rapidly 
caught  up  to  the  first  shock  wave  and  then  proceeded  as  an 
ordinary  detonation  wave  into  the  unreacted  explosive.  Separate 
equations  of  state  were  used  for  the  unrcacted  explosive  and  the 
product  gases  and,  by  assuming  additive  energies  and  volumes, 
the  equations  of  state  of  the  Intermediate  products  in  the 
reaction  zone  were  obtained.  By  stopping  the  piston  after  a 
finite  time  it  was  shown  that,  for  the  constants  used,  if  as 
little  as  one  or  two  percent  of  an  element  of  explosive  had 
reacted,  a  detonation  wave  would  subsequently  be  established, 
would  propagate,  and  could  not  be  quenched  by  any  reasonably 
strong  rarefaction-wave  which  might  appear  on  the  scene.  As 
Hubbard  and  Johnson  also  found  foi-  very  temperature  sensitive 
reactions,  the  time  necessary  for  the  major  portion  of  an 
element  of  explosive  to  react  was  small  compared  to  the  hydrody¬ 
namic  time  necessary  to  influence  motion.  However,  these 
conclusions  were  not  always  true  for  the  more  sophisticated 
numerical  experiments  to  be  described  here. 
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In  this  paper  numerical  experiments  are  described  which 
show  that  the  initiation  of  detonation  in  homoceneous  liquid  TNT 
is  a  result  of  shock  heating  and  that  the  initiation  process  is 
basically  one  of  thermal  explosion.  These  one-diir.ensional 
calculations  have  been  made  more  realistic  by  Impressing  upon 
the  liquid  TNT  acceptor  a  pressure  pulse  which  is  the  result  of 
a  Taylor  wave^  from  an  explosive  donor,  the  shock  from  which 
passes  through  an  inert  attenuator  before  reaching  the  attenua- 
tor-TNT  Interface,  The  set-up  is  shown  in  Fig.  1  where  X  =  0  is 
the  point  at  which  the  donor  explosive  is  initiated.  In  these 
calculations,  both  the  ler^gths  of  donor  explosive,  and  of 

the  gap  material,  d^,  have  been  varied  in  such  a  way  as  to  show 
the  coupled  effects  of  the  amplitude  of  the  transmitted  shock 
pressure  and  the  slope  of  the  pressure-distance  curve  behind 

6  8 

this  shock  in  the  liquid  TNT.  Previous  theoretical  calculations^^ 
assumed  that  the  state  behind  the  shock  front  was  uniform.  For 
fixed  a  critical  gap  length,  d^,  has  been  determined 

such  that  if  this  length  is  slightly  exceeded  detonation  will  not 
occur  in  the  liquid  TOT.  These  calculations  may  be  considered  as 
one-dimensional  idealized  versions  of  the  experimental  "gap  tests''^; 
the  critical  gap  length  and  resulting  transmitted  shock  amplitude 
in  the  TNT  is  then  equivalent  to  the  "505^  gap  length  and  pressure" 
values.  The  numerical  experiments  show  that  the  reaction  rate 
first  becomes  Important  at  the  attenuator-liquid  TOT  interface, 
where  the  temperature  has  been  highest  for  the  longest  time,  and 
that  detonation  first  occurs  here  even  tho’jgh  the  transmitted 
shock  may  already  have  penetrated  the  TNT  to  a  depth  of  several 
centimeters.  The  most  recent,  very  precise  physical  experiments^*^ 
on  homogeneous  explosives  support  this  conclusion. 

Section  II  contains  the  appropriate  hydrodynamic  equations 
to  be  solved.  The  equations  of  state  for  the  unreacted  explo¬ 
sive,  product  gases,  reactive  mixture,  and  gap  material  are  given 
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In  Sec.  III.  The  finite  difference  scheme  for  solving  the 
partial  differential  equations  of  hydrodynamics  is  described  in 
Sec.  IV.  Section  V  gives  the  results  of  the  one-dimensional  gap 
test  calculations  made  for  liquid  TNT.  Finally,  these  theoreti¬ 
cal  results  are  compared  to  actual  experimental  results  in 
Sec.  VI. 


II.  HYDRODYNAMIC  EQUATIONS 

The  unsteady  one -dimensional  hydrodynamic  equations  are 
most  conveniently  given  in  the  Lagrangean  form  since  this 
reveals  where  each  particle  of  fluid  came  from  initially  and 
hence,  supplies  more  information  than  the  Eulerian  form. 
Importantly,  shocks  and  moving  contact  discontinuities  can  be 
more  easily  followed,  and  mass  is  automatically  conserved.  Let 
the  Lagrangean  coordinate  x  denote  the  Eulerian  space  coordinate 
X  of  a  small  fluid  element  at  time  t  =  0,  and  X  =  X(x, t)  denote 
the  position  of  that  same  fluid  element  at  time  t  >  o,  where  x 
and  t  are  the  Independent  variables.  The  equations  representing 
the  conservation  of  mass,  momentum,  and  energy  for  a  single 
viscous  fluid  are  respectively. 


3V/3t  =  V°(x)9u/3x,  (1) 

3u/3t  =  V°{x)9(G-P)/ax  (2) 

SE/Bt  =  V°(x)(G-P)au/9x  (3) 

where 

G  =  ^[V°(x)A]  au/9x,^^  (4) 


and  E,  P,  V,  u,  and  jx  are  respectively,  the  specific  Internal 
energy,  pressure,  specific  volume,  particle  velocity,  and 
coefficient  of  viscosity.  The  superscript  o  will  always  refer 
to  initial  conditions.  The  equation  of  state. 
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E  =  r(F,V),  (5) 

then  completes  the  hydrodynamic  description,  provided  the  fluid 
species  are  unchanged.  However,  in  the  explosive  an  irreversi¬ 
ble  first-order  chemical  reaction, 

[i]  -Csl  ^  (6) 

is  postulated,  where  the  chemical  kinetic  equation  which  governs 
the  conversion  of  vinreacted  liquid  explosive  [JlJ  to  product  gas 
[s']  13 

aw/8t  =  -z  w  exp[-E*/(ffr)|.  (7) 

Here,  w,  T,  Z,  E  ,  and  R  are  respectively,  the  mass  fraction  of 
unreacted  explosive,  tenperature,  frequency  factor,  activation 
energy,  and  gas  constant.  To  determine  the  temperature  the 
additional  equation  of  state, 

E  =  s(T,V),  (8) 

must  be  given.  Equations  (l),  (2),  (5),  (4),  (5),  (7),  and 
(8)  are  seven  equations  defining  the  seven  dependent  variables 
E,  P,  V,  T,  u,  G,  and  w  in  terms  of  x  and  t.  Finally,  the 
Eulerlan  position  X(x, t)  of  the  particle  Initially  at  x  can  be 
found  from 

=  u.  (9) 

It  should  be  noted  that  E  and  V  are  the  specific  energy  and 

volume  of  an  element  of  mass  which  may  be  a  mixture  of  liquid 
and  gas.  The  element  is  assumed  to  be  a  uniform  mixture  vrlth 
uniform  pressure  and  temperature. 
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III.  EQUATIONS  OP  STATE 

For  chemical  reactions  that  Involve  a  finite  reaction  zone. 
It  Is  necessary  to  describe  not  only  the  initial  and  final 
chemical  states  but  also  the  Intermediate  ones.  This  can  be 
done  in  a  slitple  manner  if  it  is  assumed  that  E  and  V  (of  the 
intermediate  mixture)  are  equal  to  the  sum  of  the  internal 
energies  and  volumes  of  liquid  explosive  and  gaseous  product, 
l.e.. 


V  =  w\^ 

+  (l-w)Vg, 

(10) 

II 

w 

+  (l-w)Eg, 

(11) 

where  the  subscripts  and  g  refer  to  the  liquid  TNT  and  product 
gas  respectively,  and  V^,  E^,  Vg,  and  Eg  are  respectively,  the 
specific  volume  and  energy  that  the  pure  liquid  and  gas  would 
have  if  each  existed  in  isolation  at  the  pressure  and  tempera¬ 
ture  that  the  mixture  is  assumed  to  have. 

Since  the  equation  of  state  is  so  Important  in  determining 
the  fluid  flow  and  chemical  reaction  characteristics,  it  is 
desirable  to  introduce  separate  equations  of  state  for  the  liquid 
TNT  and  product  gases  rather  than  use  a  single  equation  of  state 
for  both  reactant  and  product.  For  the  product  gas,  the  equa¬ 
tions  of  state  used  are 


(12) 


S 


(15) 


For  the  liquid  TNT,  the  E, P,V  equation  of  state  is  taken  as 

E^  =  e5  +  [(P+B)V;^  -(P®+B)v/]  /{y,  -l),  (l4) 
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where 

B  =  pf(c?)2  (15) 

Is  determined  from  the  thermodynamic  identity  defining  the 
adiabatic  soxmd  speed  c, 

c2  =  [p(9P/SE)^  -  (3P/3V)j,]  ,  (I6) 

and  the  initial  energy,  E°,  is  as  yet  undetermined.  The  quanti¬ 
ties  c°  and  7^  are  evaluated  from  experimental  data  in  the 
following  manner.  From  the  three  Rankine-Hugonlot  shock  equations, 

Ph  -  V’  (P’’) 

Eg  -  Eo  =  i  (18) 

“h  ' 

and  Eq.  (l4),  it  can  be  easily  found  that 

(2U/Uj^)-  1  =  2(c^)V(UUjj)  +  (20) 

where  the  subscript  H  refers  to  values  for  the  shocked  state,  U 

is  the  shock  velocity  in  liquid  TNT,  and  P°  =  1  bar*  has  been 

12 

neglected.  Use  of  the  experimental  data  of  Gam  for 

liquid  TNT  to  111  kilobars,  and  a  least-squares  straight  line 
fit  for  (2U/Ujj)-1  vs.  2/(UUjj)  in  Eq.  (20),  gives  for  the  square 
root  of  the  slope  and  the  Intercept  respectively,  c°  =  0.2469  cm/ 
lisec  and  =  3.I78.  With  these  constants,  the  value  of  U  com¬ 
puted  from  Eq.  (20)  for  any  experimental  value  of  u^  differs 
from  the  experimental  value  of  U  by  at  most  3.6^.  Since  Gam's 
experimental  data  had  a  typical  spread  of  ±2.5^  in  shock 
velocity,  it  is  felt  that  Eq.  (14)  represents  the  Hugoniot  data 
sufficiently  accurately  for  the  computations  to  be  described. 
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Furthermore,  along  an  isentrope  through  P°,vp,  Fq.  (l4)  reduces  to 

P  =  B  [(Vj'/V^)^-  1J/y^  ,  (21) 

which  Is  the  adiabatic  Talt  equation^^  used  to  describe  water  to 
pressures  of  80  kllobars.  Equation  (l4)  and  the  assumption  that 
the  specific  heat  at  constant  volume  c^  ^ is  a  constant  defines 
the  E,V,  T  equation  of  state, 

E  =  E,°  +  c^^^T  +  [aV^^"')4»+  {'^  -1)F;  -(P°-!-B)vf  j/(Y^-l),  (2?) 

where 

That  this  is  true  can  be  readily  verified  by  noting  that  Eqs. 

(l4)  and  (22)  satisfy  the  thermodynamic  identity, 

(SE/BV)^  =  T(9P/^E)^(9E/9T)^-  P.  (24) 

Along  an  Isentrope, 

T  =  kv^"T:,  (25) 

where  K  is  a  constant  that  depends  on  the  particular  isentrope. 

Let  Q  be  the  internal  energy  difference  between  reactant  and 
product  at  the  initial  pressure  and  temperature,  i.e., 

Q  =  E^(P°,T°)  -  Eg(P°,T°).  (26) 

Substitution  of  Eqs,  (l;)  and  (22)  into  (26)  then  defines  E°  as 

E°  =  Q  +  c^^gT®  .  (27) 
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The  inert  gap  material  is  assumed  to  have  a  P,V  Hugoniot 
curve  given  by^ 

P  =  58.4(0.847  -  Vq)/(0.272  +  ,  (28) 

where  P  and  Vq  are  in  kilobars  and  cm'/g,  respectively;  the 
isentropic  expansion  curve  is  assumed  to  follow  the  Hugoniot. 

It  was  shown  by  Walsh  and  Christian  ^  that  the  use  of  the  P,V 
Hugoniot  to  represent  the  isentropic  e:5>anslon  curve  for  any 
material  is  equivalent  to  minimizing  u^g  -  Ujj,  where 


o 


Ufg  is  the  free  surface  velocity  at  P  =  0,  and  Uj.g  -  Ujj  is  the 
increment  in  particle  velocity  due  to  a  centered  single  rare¬ 
faction  wave  from  the  shocked  state  to  0,uj.g. 

The  following  constants  are  used  for  each  fluid; 

Tetryl(donor  exploslve)^^;  P®  =  1  bar,  p°  =  I.63  g/cm^, 

T°  =  500OK,  Q  =  1047  cal/g 

Tetryl  product  gas^^;  =  2,69 

Liquid  TNr(acceptor  e:qploslve) :  P°  =  1  bar, 

V?  0,679  cm^/g^^,  =  5540K^^ 

Q  =  856  cal/g,  =  5.178,  c°  =  0.2469  cm/psec 
c^^,  =  0.528  cal/(g  deg)^^,  B  =  89.805  kbar, 

A  =  1.6998  kbar  (cm^/g)^*^^® 

E*  =  54.4  kcal/mole^®,  Z  =  10^^*^  sec"^ 

The  values  of  Q  =  856  cal/g  and  =  2.689  for  the  detona¬ 
tion  of  liquid  TNT  were  computed  by  solving  for  Q  and  y  from  the 

s 

equations. 


2(7|-1)Q  = 
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which  are  valid  for  a  Chapman -Jovjget  detonation,  vrtiere  the 

detonation  velocity,  D  =  0.66  cn\/M-sec,  was  measured  by 

L.  B.  Seely,  Jr,^^,  and  =  I7I.8  kbars  was  determined  by 
20 

Gam  .  The  initial  temperature  and  density  in  their  work 
were  respectively,  T°  =  365®K  and  p?  =  1.455  g/cm^. 

IV.  FINITE  DIFFERENCE  EQUATIONS 

In  order  to  carry  out  the  numerical  solution  of  Eqs.  (l)  - 
(5),  and  (7)  -  (9),  they  are  replaced  by  an  equivalent  system 
of  finite  difference  equations.  The  exact  form  of  these 
equations  is  governed  by  stability  considerations  and  a 
stability  analysis  as  proposed  by  von  Neumann'  has  been 
carried  out  in  order  to  verify  that  small  errors  will  not  grow 
during  the  conputations .  The  difference  equations  used  for  the 
numerical  calculations  are  as  follows: 


(50) 

-  u5:|) , 

(51) 

(32) 

^n+1  8.  „0/  n++  n-(4\ //„n+l  ,.,n\ 

(33) 

(34) 

^  -  wj-i"  -iZ4t"(w5'*i+vi^'4)e*p  [-EV(irrpJ  , 

(35) 

(36) 
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c"i 

,n+l 


v-n-i-  i»*.n/  n+4  ,  n-i\ 

x"  .  +  u"4)  , 

J  J"^  ^  '2 


At^  =  =  i(v5+i  +  vp. 


(37) 

(38) 

(39) 


where 


riH-i 


J  =0,1,2,...,  J;  n  =  0,1,2,... 


u".  f  =  u(x4-^AX,  t-4-g-At),  w^~^  =  w(x+AX,  t-|-At) ,  etc. 

J'^2  J 

Here  x  Is  defined  as  shovm  In  Fig.  1  and  ax  (v'hlch,  vflth  more 
precision  but  also  more  rotational  complication,  should  be 
written  as  ax^)1s  the  mesh  size  in  the  1th  fluid,  where 
1  =  1,  2,  5  refers  to  respectively,  the  donor  explosive,  inert 
gap,  and  acceptor  explosive  (liquid  TNT).  The  total  number  of 
mesh  points  used  In  the  computations  was  IJO,  i.e,,  J  =  129. 

If  all  the  variables  are  known  at  t’^"^  ^^d  t”^,  the  times 

\ 

corresponding  to  the  nth  time  step  (t“'2‘<o,  t°  =  O),  and  the 
outer  boundary  conditions  are  prescribed  at  ^  =  0  and  J  =  J, 
then  the  values  at  the  new  times,  t’^’^  =  and  t^'*’^  = 

t’^+At^'*'^  are  computed  as  follows.  Let  ^  i,  J  =  '^■^>'^2* 

0<I^<l2<J,  instantaneous  positions  of  the  donor  explosive- 

inert  and  inert  -acceptor  explosive  interfaces  (contact  discon¬ 
tinuities)  respectively.  Therefore  Uj  ^  1,  J  = 
velocities  of  these  Interfaces.  These  velocities  must  be 
conputed  from  the  special  formula. 


n+i  n— ^  ,  /.4.n/,„\  ^n 

\4  =  +  ('t  /«)  5  i+i-  I  =  h  4 


where 


N 

/ 


-  v°(x)a(G-p)/^x  . 


(40) 


(41) 
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FIG.  I  IDEALIZED  ONE-DIMENSIONAL  GAP  TEST 
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The  unknovm  ^  Is  obtained  by  expanding  ^  to  first  order 
about  the  points  at  I,  I+l,  1+^  and  solving  the  resulting 

four  equations  In  the  four  \inlcnovms, 

^  ^I+i* 

Uj^,  0<J^J-1,  then  computed.  Equations  (35)^  (51), 

and  (33)  then  give  and  I^J^J-l. 

The  time  step  At”"^  Is  calculated  from 


At""*^  -  |^(Atj'^)j^jjj,^j  ,(AtJ)^j^^jij  J  .  (42) 

At  every  point  (including  those  In  the  reaction  zone  where 
0<Wj"^<l)  the  hydrodynamic  time  step  (At”''^)hy(lj.o  1b  determined 
from  the  stability  criterion. 


i(v5)^(a5-K}J“i) 


However,  in  the  reaction  zone,  o”  emd  |^(ap/aE)^*'  are  coinputed  by 

asBumlng  that  w  ->  1.  This  is  nmch  slirpler  than  talcing  Into 
account  the  change  In  w  and  has  worked  very  well  without  any  loss 
In  computational  time.  In  the  reaction  zone  the  critical  chemical 
time  step  is 


(At”)  -(2/Z)exp  [eV(Rt5)]  ,  (44) 

"  chem ■  " 

which  is  derived  from  Eq,  (35)  by  letting  Wj"^^  ■  0. 

For  unreacted  liquid  TEfl  and  product  gas,  Eq.  (34)  is  defined 
by  Eqs,  (12)  and  (l4)  in  which  e”"^^  is  linear  in  Hence, 

Eqs.  (32)  and  (34)  can  bo  solved  explicitly  for  and 
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The  teinperature  then  follows  from  Eq.  (56)  (See  Eqs,  (I5) 

and  (22)).  For  the  gap  material  is  computed  from  Eq.  (28), 

and  follows  from  Eq.  (32) .  However,  for  the  reactive  mix¬ 

ture  the  procedure  is  much  more  complex.  Equations  (lO),  (11), 
(12),  (13),  (l4),  and  (22)  evaluated  at  J,  n+1  and  Eq,  (32)  are 
seven  equations  in  the  nine  variables 

and  w^"^.  But  and  w”+4  are  already 

known  so  that  the  first  seven  variables  listed  can  be  deter¬ 
mined.  However,  this  first  Involves  the  extraction  of  the  root 
of  the  equation 


^(X)  =  (bj^+b2X+b^X^)x'^V  +  (bi^+b^X)X  =0,  (45) 

which  is  carried  out  to  five  place  accuracy.  Here,  b^,  bg,  b^, 
bi^,  and  b^  are  functions  of  the  flow  variables  which  are 
already  known  at  n  or  n+1. 

During  the  progress  of  every  machine  computation,  the 
conservation  of  mass,  momentum,  and  energy  were  continually 
checked.  Mass  and  momentum  were  extremely  well  conserved; 
energy  losses  were  at  most  1  or  2^, 

V.  ONE -DIMENSIONAL  GAP  TEST  CALCULATIONS 
The  one -dimensional  idealized  version  of  the  gap  test  is 
shown  in  Fig.  1.  A  nximber  of  numerical  e:q>erlraents  were 
performed  on  an  IBM  7090  computer  in  which  the  length  of  the 
donor  e:q)losive  (Tetryl),  was  fixed  and  the  length  of 

the  gap  d^  was  varied  until  detonation  in  the  liquid  TNT  was 
Just  barely  possible.  For  gaps  larger  than  this  critical 
length  dg,  detonation  never  propagated;  for  smaller  gaps 
detonation  always  propagated.  The  results  of  these  machine 
computations  are  given  in  Table  1.  Here  P^  is  the  value  of 
the  pressure  of  the  transmitted  shock  as  it  enters  the  liquid 
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TMT.  For  Cases  1  -  5,  the  shock  initiated  reaction  upon  enter¬ 
ing  the  liquid  TUT  and  had  advanced  a  distance  x  into  the 
explosive  at  the  time  the  explosive  cell  at  the  gap-explosive 
interface  detonated.  For  Cases  2,  5,  and  4,  the  shock  anplltude 
continuously  decreased  for  t<'T  •  The  induction  time  t  is  the 
"cooking"  time  of  that  cell  after  the  shock  passed  over  it. 

The  pressure  distribution  at  various  times  is  shovm  in 
Fig.  2  for  Case  2.  In  Fig.  2(a)  the  Taylor  detonation  wave  has 
already  reflected  off  the  product  gas-attenuator  interface  and 
a  reflected  rarefaction  vfave  is  moving  back  into  the  product 
gas  while  a  shock  is  transmitted  into  the  gap  material.  Only 
the  rightmost  part  of  the  product  gas  is  shown  since  most  of 
the  gas  is  expanding  to  the  left  into  a  vacuum  (approximating 
air  at  several  bars).  In  Fig.  2(b)  the  progress  of  the  shock 
in  the  gap  is  sho\m;  and  in  Fig.  2(c)  the  shock  has  reflected 
off  the  gap-liquid  TIJT  interface  with  the  result  that  higher 
an^lltude  shocks  are  moving  into  the  gap  and  into  the  explosive. 
In  Fig.  2(d)  the  almost  discontinuous  rise  in  pressure  at  the 
interface,  due  to  reaction,  can  be  seen.  Now,  in  a  very  short 
time  coirplete  reaction  occurs  as  shown  in  Fig.  2(e).  The 
detonation  wave,  traveling  faster  than  the  forward  shock  front, 
rapidly  catches  up  to  the  latter  (Figs.  2(f)  and  2(g)),  and 
then  proceeds  as  an  ordinary  detonation  wave  (Pig.  2(h)).  In 
Figs.  2(d)  to  2(g),  it  can  be  seen  that  the  amplitude  of  the 
transmitted  shock  is  decreasing  because  of  rarefaction  waves 
from  the  rear.  Though  the  pressure  at  the  gap-liquid  TNT 
interface  is  lower  than  at  the  shock  front  the  temperature  is 
always  higher  because  of  the  release  of  energy  in  the  conver¬ 
sion  of  explosive  to  product  gas.  In  Fig.  3  the  pressures  at 
the  interface  and  at  the  shock  front  are  plotted  as  functions 
of  time  for  Case  2, where  the  time  scale  has  been  shifted  such 
that  the  new  origin  coincides  with  that  (real)  time  at  which 
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Pig.  2  The  pressure  distribution  In  the  one- 

dimensional  gap  test  for  various  values 
of  t  for  Case  2.  The  heavy  vertical  lines 
A,  B,  C,  and  D  represent  respectively,  the 
air-product  gas  interface,  the  product  gas- 
attenuator  Interface,  the  attenuator-liquid 
TNT  (or  liquid  TNT  product  gas)  interface, 
and  the  liquid  TNT-air  Interface,  Each  dot 
represents  a  mass  point  whose  initial 
position  is  given  by  the  Lagrangean  coordi¬ 
nate  X.  The  motion  of  any  mass  point  in 
the  X  direction  can  be  followed  by  observing 
the  motion  of  that  same  dot  (provided  the 
observer  has  a  pair  of  telescopic  eyes). 

The  values  of  t  (nsec),  as  measured  from 
the  instant  at  which  the  shock  first  entered 
the  explosive,  are  as  follows:  (a)  -7.95» 
(b)  -2.80,  (c)  0.45,  (d)  5.84,  (e)  4.6?, 

(f)  6.15,  (g)  7.85,  (h)  10.74. 
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the  shock  first  entered  the  explosive.  Complete  reaction  occurs 
at  T  =  4.52  lisec.  In  Case  2  the  rarefaction  wave  attenuates  the 
shock  that  has  progressed  Into  the  liquid  TNT  Interior,  and 
chemical  reaction  at  this  front  Is  never  able  to  compete  success¬ 
fully  with  the  rarefaction  wave.  This  wave  effectively  annihi¬ 
lates  any  compression  waves  which  may  originate  at  the  slowly 
reacting  Interface  and  which  would  tend  to  move  In  the  direction 
of,  and  thus  reinforce,  the  shock  front. 

For  Cases  1  and  5>  the  results  are  similar  to  Case  2, 
except  that  at  the  gap-llquld  TNT  Interface  the  pressure  (and 
temperature)  continuously  Increases  with  time  during  reaction 
as  shown  In  Fig.  5  for  Case  5.  Therefore,  the  pressure  and 
temperature  decrease  monotonlcally  In  the  direction  of  the  shock 
front.  In  these  cases  It  Is  not  possible  to  decide  whether  the 
pressure  at  the  shock  front  decreases  or  Increases  with  time  for 
t  <  T  since  the  numerical  method  does  not  define  a  sharp  shock 
front  and,  hence.  It  Is  not  possible  to  evaluate  accurately  the 
pressure. 

It  Is  at  the  Interface  that  the  Interaction  between  the 
shock  pressure  amplitude  and  the  varying  pressure  distribution 
behind  the  shock  Is  seen.  In  previous  work,  the  state  behind 
the  Initial  shock  was  uniform,  thus  leading  to  an  increase  In 
both  P  and  T  at  the  piston-explosive  Interface^  or  at  the  Inert 

Q 

gap-explosive  interface  .  However,  In  Case  2,  the  pressure  at 
the  Interface  decreases  with  time  during  most  of  the  "cooking" 
period  because  of  the  rarefaction  wave  from  the  rear.  Increas¬ 
ing  the  transmitted  shock  pressure  and  the  pressure  gradient  at 
the  interface,  as  In  Case  1,  now  leads  to  a  monotonlc  pressure 
Increase  with  time  at  the  interface.  Thus,  here  the  effect  of 
pressure  anplltude  is  more  Important  in  increasing  the  reaction 
rate  (through  higher  tenperatures)  than  that  of  Increased 
pressure  gradient  in  decreasing  the  reaction  rate  (throxigh 
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FIG.  3  PRESSURES  AT  THE  ATTENUATOR-LIQUID  TNT  INTERFACE  AND  AT  THE  SHOCK 
FRONT  FOR  CASES  2  AND  5 


NOLTR  62-159 


lowering  the  temperature) . 

The  t,X  diagram  is  particularly  informative  for  showing 
the  wave  and  particle  motion  and  this  is  depicted  in  Pig,  4  for 
Case  1.  Only  motion  in  the  vicinity  of  the  gap-liquid  TNT 
interface  is  shown.  Curve  C^  represents  the  motion  of  the  shock 
wave  in  the  gap  before  it  strikes  the  gap-explosive  interface, 

Cg  represents  the  path  of  the  interface  which  moves  to  the  right 
with  the  local  particle  velocity  while  a  transmitted  shock  wave 
with  decreasing  velocity  moves  along  C^  into  the  explosive. 
Actually,  for  the  first  part  of  its  travel  along  Cy  the  shock 
front  immediately  initiates  chemical  reaction.  High-order 
detonation  originates  at  the  Interface  after  a  time  t  and  pro¬ 
ceeds  through  the  pre-shocked  and  very  slightly  reacted  material 
(0  to  5  percent  reacted)  along  C^  until  the  detonation  overtakes 
the  transmitted  chock  at  the  transition  point  A  and  passes  into 
the  unshccked  explosive  along  C^.  A  similar  description  holds 
for  Cases  1  and  5  (except  perhaps  for  the  shock  velocity  in  the 
liquid  TNT  before  detonation). 

Interesting  results  were  obtained  in  Case  5  with  d^^  =  2.95 
cm  and  v:lth  =  2.54  cm  (i.e,,  a  perturbation  of  Case  2). 

As  in  Case  2,  the  first  explosive  cell  detonated,  but  subse¬ 
quently  the  detonation  was  quenched  by  the  rarefaction  v^ave 
from  the  rear  and  no  more  than  that  one  cell  ever  reacted  to 
completion.  This  sequence  of  events  is  shown  in  Fig.  5.  A 
slight  change  to  d^  =  2.94  cm  (i.e..  Case  4)  resulted  in  failure 
to  detonate  at  any  point.  In  other  words,  d^  =  2.92  cm  can  be 
very  accurately  determined. 

VI.  COMPARISON  OF  THEORETICAL  WITH  EXPERIMENTAL  RESULTS 

Extensive  experiments  on  the  shock  initiation  of  detonation 
in  homogeneous  explosives  were  performed  by  Campbell,  Davis,  and 
Travls^^.  Most  of  the  work  was  done  on  nltromethane,  but  liquid 
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FIG.  4  t,  X  DIAGRAM  SHOWING  THE  WAVE  AND  PARTICLE 
MOTION  FOR  CASE  2 
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Fig.  5  The  pressure  distribution  in  the  one¬ 

dimensional  gap  test  for  various  values 
of  t  for  Case  3.  Reaction  begins  at 
the  attenuator-liquid  TNT  Interface  In 

(a) j  complete  reaction  at  Interface  In 

(b) ;  quenching  of  the  reaction  by  rare¬ 
faction  wave  In  (c)  and  (d).  See  Pig.  2 
for  definition  of  the  lines  A,  B,  C,  and 
D.  The  values  of  t  (p-sec),  as  measured 
from  the  Instant  at  which  the  shock 
first  entered  the  explosive,  are  as 
follows;  (a)  4.45,  (b)  5.40,  (c)  6.36, 

(d)  10.95. 
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TNT,  liquid  DINA,  Dlthekete  15,  and  single  crystals  of  PETN  were 
shown  to  behave  in  the  same  way.  Plane  wave  lenses  were  com¬ 
bined  with  additional  high  explosive  and  an  inert  shock -attenua¬ 
tor  to  produce  a  shock  wave  of  desired  anplitude  in  the  acceptor 
explosive.  The  resulting  t,X  representation  of  the  Initiation 
behavior  of  the  homogeneous  explosives  is  qualitatively  similar 
to  that  shown  in  Fig.  4  .  Their  experimental  results  show  that 

the  initial  shock  wave  in  a  homogeneous  explosive  has  a  constant 
or  slightly  decaying  velocity  as  a  function  of  time;  the  compu¬ 
tations  for  Case  2  show  a  slightly  decaying  shock  velocity  as 
evidenced  by  the  curvature  of  in  Fig,  4  or  by  the  decreasing 
shock  pressure  in  Fig.  5.  The  experimental  results  show  that 
small  changes  in  the  initial  transmitted  shock  pressure  P^  lead 
to  large  changes  in  the  induction  time  t ;  the  theoretical  conpu- 
tatlons  listed  in  Table  I  show  the  same  characteristic.  The 
experimental  work  indicates  that  the  detonation  wave  overtakes 
the  transmitted  shock  and  tenporarlly  overdrives  detonation  in 
the  unshocked  explosive  ahead  of  it.  For  nitromethane  the 
detonation  velocity  increase  is  approximately  10^  and  the  veloc¬ 
ity  decreases  in  a  few  microseconds  to  normal  detonation  velocity; 
the  calculations  show  a  slight  overshoot  along  in  Fig.  4  that 
disappears  after  about  1  p,sec.  The  preliminary  initiation 
data^*^  for  liquid  TNT,  initially  at  558°K,  gave  as  the  initiating 
pressure,  125  kllobars,  and  as  the  induction  time,  0.7  p,sec. 

This  compares  with  a  pressure  of  8l  -  89  kllobars  and  Induction 
times  of  several  ^.secs  from  Table  I.  The  discrepancy  in  the 
initiating  pressures  between  the  e^qjerlmental  and  theoretical 
results  may  be  due  to  the  equation  of  state,  Eq,  (l4),  which 
represents  the  experimental  U,Ujj  Hugoniot  relatively  faithfully, 
but  which  can  lead  to  large  errors  in  the  temperature  (and 
throvigh  the  chemical  kinetics,  to  large  errors  in  the  Initiating 
pressure)  if  is  not  really  a  constant  as  assumed  in  the 
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TABLE  I.  Summary  of  calculational  results  for  liquid  TNT,  x_ 
is  the  distance  from  the  initial  position  of  the  gap -liquid  TNT 
Interface  to  the  shock  front  position  when  complete  reaction 
first  occurred  at  the  Interface  (see  Pig.  4).  t  is  the  time 
Interval  between  the  time  that  the  shock  first  entered  the 
liquid  TNT  and  the  time  complete  reaction  first  occurred  at  the 
gap-liquid  TNT  interface  (see  Pig.  4). 


Case 

‘^donor 

cm 

cm 

kbar 

^s 

cm 

T 

p,sec 

1 

2.54 

2.70 

85 

0.7 

1.5  ® 

2 

2.54 

2.92  ^ 

81 

1.9 

4.5 

5 

2.54 

2.93  ^ 

80<P<8l 

2.6 

6.4 

4 

2.54 

2.94  ° 

80<P<8l 

00 

00 

5 

15.24 

15.1  ^ 

89 

1.5 

2.2  ® 

The  first  liquid  TNT  mass  cell  went  to  complete  reaction  but 
the  detonation  was  subsequently  quenched  by  the  rarefaction 
wave  and  no  other  explosive  cells  went  to  complete  reaction. 

°  19^  of  the  first  liquid  TNT  mass  cell  reacted.  The  rarefac¬ 
tion  wave  subsequently  quenched  all  further  reaction  and 
detonation  was  never  initiated  in  Case  4. 

°  dg  <  d^  for  this  Case. 

®  Why  is  T  greater  for  Case  5  than  for  Case  1  since  is 

greater  and  the  pressure  gradient  is  smaller  at  the  gap -liquid 
TNT  Interface?  This  is  probably  due  to  the  difference  in  the 
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TABLE  I  (Continued) 

viscosities  p.  of  the  gap  material  between  Cases  1-4  (the 
same)  and  Case  5.  While  the  numerical  scheme  spreads  the 
major  part  of  the  shock  front  over  5  or  6  cells,  there 
always  exists  a  low  pressure  perturbation  that  moves  ahead 
of  the  main  shock  front  and  hence  slightly  changes  the 
initial  conditions  that  the  shock  front  sees.  These  pertur¬ 
bations  are  of  slightly  different  magnitude  for  Cases  1  and 
5.  The  result  is  that  the  shock  temperature  does  not  con¬ 
verge  as  quickly  to  its  true  value  (an  upper  bound  for  the 
numerical  temperature)  for  Case  5  as  it  does  for  Case  1. 

This  fact  explains  the  Inversion  of  the  values  of  ,  v;hich 
is  a  very  sensitive  function  of  temperature. 

It  should  be  mentioned  that  the  precursor  shock,  caused 
by  use  of  an  artificial  viscosity,  is  minimized  by  testing 
the  relative  change  in  the  volume  at  each  time  step  and  by 
bypassing  the  hydrodsmamlc  calculations  if  this  relative 
volume  change  is  less  than  10"  . 
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derivation  of  Eq.  (22)  .  Also,  errors  in  the  evaluation  of 

E  and  Z,  or  changes  in  these  values  with  Increased  P  and  T,  can 
easily  force  the  initiating  temperature  and  (through  the  equa¬ 
tion  of  state)  the  initiating  pressure  upwards.  Furthermore, 
the  pressure  distribution  behind  the  initiating  shock  in  the 
liquid  TNT  is  surely  not  the  same  in  the  experimental  work  and 
in  the  theoretical  calculations.  This  may  account  for  some  of 
the  difference.  Since  the  calculated  induction  time  is  of  the 
correct  order  of  magnitude,,  and  particularly  since  the  t,X  wave 
diagram  for  Case  2  is  in  qualitative  agreement  with  experimental 
observations,  it  is  felt  that  the  model  used  in  the  above  calcu¬ 
lations  is  reasonable. 
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